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DE TECTION, IM AGINCLAN Q CHARACTERIZA TTONjjFRTRF. A ST TUMORS 



This application claims priority from U.S. Provisional Application Serial No. 
60/074,504 filed on February 1 1, 1998 and from U.S. Provisional Application Serial No. 
5 60/098,01 8 filed on August 26, 1998, both of which are incorporated by reference as if fully 
set forth herein. 

The Field of the Invention 

The present invention relates to imaging and qualitative or quantitative 
10 characterization of biological tissue using visible or infra-red radiation, and more particularly 
to detection, imaging and characterization of breast tumors. 

Background 

Traditionally, potentially harmful ionizing radiation (for example X-ray or y-ray) has 
15 been used to image biological tissue. This radiation propagates in the tissue on straight, 
ballistic tracks, i.e., scattering of the radiation is negligible. Thus, imaging is based on 
evaluation of the absorption levels of different tissue types. For example, in roentgenography 
the X-ray film contains darker and lighter spots. In more complicated systems, such as 
computerized tomography (CT), a cross-sectional picture of human organs is created by 
20 transmitting X-ray radiation through a section of the human body at different angles and by 
electronically detecting the variation in X-ray transmission. The detected intensity 
information is digitally stored in a computer which reconstructs the X-ray absorption of the 
tissue at a multiplicity of points located in one cross-sectional plane. 

Near infra-red radiation (NIR) has been used to study non-invasively the oxygen 
25 metabolism in tissue (for example, the brain, finger, or ear lobe). Using visible, NIR and 

infra-red (IR) radiation for medical imaging could bring several advantages. In the NIR or IR 
range the contrast factor between a tumor and a tissue is much larger than in the X-ray range. 
In addition, the visible to IR radiation is preferred over the X-ray radiation since it is non- 
ionizing and thus, potentially causes fewer side effects. However, the visible or IR radiation 
30 is strongly scattered and absorbed in biological tissue, and the migration path cannot be 
approximated by a straight line, making inapplicable certain aspects of cross-sectional 
imaging techniques. 



WO 99/40840 PCT7US99/02953 

-2- 

Computerized Tomography using NIR spectrometry has been used for in vivo 
imaging. This technique utilizes NIR radiation in an analogous way to the use of X-ray 
radiation in an X-ray CT. The X-ray source is replaced by several laser diodes emitting light 
in the NIR range. The NIR-CT uses a set of photodetectors that detect the light of the laser 
diodes transmitted through the imaged tissue. The detected data are manipulated by a 
computer similarly as the detected X-ray data would be in an X-ray CT. Different NIR-CT 
systems have recognized the scattering aspect of the non-ionizing radiation and have modified 
the X-ray CT algorithms accordingly. 

The above-mentioned X-ray or y-ray techniques have been used to detect a tissue 
tumor. Under the term "angiogenesis" I mean the generation of new blood vessels into a 
tissue or organ. Under normal physiological conditions humans or animals undergo 
angiogenesis only in very specific restricted situations. For example, angiogenesis is 
normally observed in wound healing, fetal and embryonal development and formation of the 
corpus luteum, endometrium and placenta. 

Both controlled and uncontrolled angiogenesis are thought to proceed in a similar 
manner. Persistent, unregulated angiogenesis occurs in a multiplicity of disease states, tumor 
metastasis and abnormal growth by endothelial cells and supports the pathological damage 
seen in these conditions. The diverse pathological disease states in which unregulated 
angiogenesis is present have been grouped together as angiogenic dependent or angiogenic 
associated diseases. The hypothesis that tumor growth is angiogenesis dependent was first 
proposed in 1971. (Folkman J., Tumor angiogenesis: Therapeutic implications., N. Engl. 
Jour. Med. 285: 1 182-1186, 1971) In its simplest terms it states: "Once tumor 'take' has 
occurred, every increase in tumor cell population must be preceded by an increase in new 
capillaries converging on the tumor." Tumor 'take' is understood to indicate a prevascular 
phase of tumor growth in which a population of tumor cells occupying a few cubic 
millimeters volume and not exceeding a few million cells, can survive on existing host 
micro vessels. Expansion of tumor volume beyond this phase requires the induction of new 
capillary blood vessels. This explanation was directly or indirectly observed and documented 
in numerous publications. 
) Breast cancer is among the most common and the most feared malignancies in women. 

It has an unpredictable course, the treatment is frequently physically and emotionally draining 
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and the risk of metastatic spread persists for many years. Due to its high occurrence rate, 
routine breast cancer screening, which includes physical examination and x-ray 
mammography, plays an important role in current health care. X-ray mammography can 
detect perhaps 90% of all masses and increases the 10-year survival rate to about 95% for 
patients with cancers solely detected by mammography. Although the modern mammography 
uses a low-dose of x-rays, it still involves some very small risk of inducing cancers by the 
radiation. Other tests, such as magnetic resonance imaging (MRI) and gadolinium enhanced 
MRI, have been used successfully for detection of breast tumors and may be used routinely 
for screening in the future. 

After a small suspicious mass is detected in the breast non-invasively, an excisional 
biopsy is usually performed to exclude or diagnose malignancy. The biopsy specimen is 
removed under local anesthesia and is used for histopathological diagnosis. The statistics 
show that up to 75%of the excisional biopsies, the biopsied tissue is diagnosed to be benign. 
Thus, a majority of patients undergoes this unpleasant and costly procedure unnecessarily. 

Therefore, a non-invasive, relatively inexpensive technique that can detect and 
characterize breast tumors may find its place in today's health care alone or in conjunction 
with the above-mentioned techniques. 



Summary 

The present invention includes different novel apparatuses and methods for optical 
examination of biological tissue, in general, and breast tissue, specifically, using visible or 
infra-red light. The optical examination technique can be used alone to detect and 
characterize a breast tumor, or can be used in combination with X-ray mammography, 
ultrasound examination,yMRI, or a needle biopsy. Furthermore, the optical examination 
technique can be used to examine women of any age. 

The technique can employ one or several optical modules positioned on the right or 
left breast of a female patient placed in different positions. The patient may be sitting upright 
supporting the examined breast by the optical module, or may be lying face down with the 
breast on the optical module pad. Alternatively, the patient may be lying supine face up with 
the breast spread over the chest as evenly as possible. If a suspicious mass is detected, the 
technique can non-invasively characterize the mass by taking optical data at different 



WO 99/40840 PCT/US99/02953 



wavelengths and by measuring one or several tissue specific characteristics related to the 
tissue metabolism (or hypermetabolism), biochemistry, pathophysiology (including 
angiogenesis) or another characteristic of a pathological tissue condition. 

In one aspect, the optical examination technique employs an optical system for in vivo 

5 non- invasive examination of a volume of biological tissue of a subject. The optical 

examination system includes an optical module, a controller and a processor. The optical 
module includes an array of optical input ports and optical detection ports located in a selected 
geometrical pattern to provide a multiplicity of source-detector paths of photon migration 
inside the biological tissue. Each optical input port is constructed to introduce into the tissue 

10 volume visible or infrared light emitted from a light source. Each optical detection port is 
constructed to provide light from the tissue to a light detector. The controller is constructed 
and arranged to activate one or several light sources and light detectors so that the light 
detector detects light that has migrated over at least one of the source-detector migration 
paths. The processor receives signals corresponding to the detected light and creates a defined 

1 5 spatial image of the examined tissue. 

The optical examination system may generate single wavelength or multiple 
wavelength images of the examined tissue, wherein the employed wavelength is sensitive to 
absorption or scattering by a tissue constituent (e.g., an endogenous or exogenous pigment, 
tissue cells, chemical compounds) or is sensitive to structural changes in the tissue. The 

20 optical images may display tissue absorption, tissue scattering, or both. The optical imaging 
system may also generate blood volume hemoglobin oxygenation images and hemoglobin 
deoxygenation images (or images of any other tissue constituent) based on a single 
wavelength optical data or a multiple wavelength optical data. A processor may use different 
image processing and enhancing algorithms known in the art. The processor may correlate 

25 several images taken on the same tissue or taken on symmetrical tissue regions such as the left 
breast and the right breast, or the left arm and the right arm. Based on this correlation, the 
system detects a suspicious tissue mass and characterize the detected mass. The correlation 
includes determining congruency of the structures detected in different images. The processor 
may employ different types of combined scoring, based on several optical images alone or in 

30 combination with X-ray mammography, ultrasound examination, or ./MRI, to characterize a 
suspicious tissue mass. 
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The optical imaging system may generate the above-described images by examining 
symmetrical tissue regions of the right breast and the left breast, or may generate images of 
both the entire right breast and the entire left breast. To identify and characterize a suspicious 
tissue mass, the processor may employ the different types of combined scoring by correlating 
5 the right breast image with the left breast image. 

The optical imaging system may collect single wavelength or multiple wavelength 
data of a breast tissue model for calibration or for detection of background data. In the 
calibration procedure, the optical module is placed on the model, and the imaging system can 
collect a limited number of optical data or can collect optical data using the same sequences as 
10 used during the tissue examination. The system may either collect and store the model data 
for subsequent digital processing, or may adjust the source or detector gains to detect optical 
data according to a predetermined optical pattern. The imaging system may use different 
breast models having the same scattering coefficient and the same absorption coefficient as 
the normal breast tissue of a female, for example, below 40 years or above 40 years. 
15 Furthermore, the models may have different sizes and have a shape of a female breast during 
examination. For example, the model has a shape of the breast (or the breast and a portion of 
the chest) of a woman lying on her back during examination. 

According to another aspect, the optical examination technique employs an optical 
system for in vivo, non-invasive examination of biological tissue of a subject. The optical 
20 system includes an optical module, a controller , and a processor. The optical module 

includes an array of optical input ports and detection ports located in a selected geometrical 
pattern to provide a multiplicity of photon migration paths inside an examined region of the 
biological tissue. Each optical input port is constructed to introduce visible or infrared light 
emitted from a light source. Each optical detection port is constructed to receive photons of 
25 light that have migrated in the examined tissue region from at least one of the input ports and 
provide the received light to a light detector. The controller is constructed and arranged to 
control operation of the light source and the light detector to detect light that has migrated 
over at least one of the photon migration paths. The processor is connected to receive signals 
from the detector and arranged to form at least two data sets, a first of the data sets 
30 representing blood volume in the examined tissue region and a second of the data sets 

representing blood oxygenation in the examined tissue region. The processor is arranged to 



W ° 99/40840 PCT/US99/02953 

-6- 

correlate the first and second data sets to detect abnormal tissue in the examined tissue region. 

Preferably, the second data set includes hemoglobin deoxygenation values. The 
processor may be arranged to form a third data set being collected by irradiating a reference 
tissue region. 

According to another aspect, the optical examination technique employs an optical 
system for in vivo, non-invasive examination of biological tissue of a subject. The optical 
system includes an optical module, a controller , and a processor. The optical module 
includes an array of optical input ports and detection ports located in a selected geometrical 
pattern to provide a multiplicity of photon migration paths inside an examined region of the 
biological tissue. Each optical input port is constructed to introduce visible or infrared light 
emitted from a light source. Each optical detection port is constructed to receive photons of 
light that have migrated in the tissue from at least one of the input ports and provide the 
received light to a light detector. The controller is constructed and arranged to control 
operation of the light source and the light detector to detect light that has migrated over at 
least one of the photon migration paths. The processor is connected to receive signals from 
the detector and arranged to form at least two data sets, a first of the data sets being collected 
by irradiating an examined tissue region of interest and a second of the data sets being 
collected by irradiating a reference tissue region having similar light scattering and absorptive 
properties as the examined tissue region. The processor is arranged to correlate the first and 
second data sets to detect abnormal tissue in the examined tissue region. 

According to another aspect, the optical examination technique employs an optical 
system for in vivo, non-invasive examination of biological tissue of a subject. The optical 
system includes an optical module, a controller , and a processor. The optical module 
includes an array of optical input ports and detection ports located in a selected geometrical 
pattern to provide a multiplicity of photon migration paths inside an examined region of the 
biological tissue or a model representing biological tissue. Each optical input port is 
constructed to introduce visible or infrared light emitted from a light source. Each the optical 
detection port is constructed to receive photons of light that have migrated in the tissue or the 
model from at least one of the input ports and provide the received light to a light detector. 
The controller is constructed and arranged to control operation of the light source and the light 
detector to detect light that has migrated over at least one of the photon migration paths. The 



WO 99/40840 



PCT/US99/02953 



-7 - 



processor is connected to receive signals from the detector and arranged to form at least two 
data sets of two tissue regions, a first of the data sets being collected by irradiating an 
examined tissue region and a second of the data sets being collected by irradiating a region of 
a tissue model having selected light scattering and absorptive properties. The processor is 
5 arranged to correlate the first and second data sets to detect abnormal tissue in the examined 
tissue region. 

Preferred embodiments of these aspects of the inventions have one or more of the 
following features. 



10 determining congruence between data of the two sets. 

The processor may be programmed to order the first and second data sets as two- 
dimensional images and to determine the congruence using the two-dimensional images. 
The processor may be programmed to order the first and second data sets as two-dimensional 
images and to determine the congruence using the following formula: 



The processor may be further arranged to determine a location of the abnormal 
tissue within the examined tissue region. 

The processor may be adapted to produce from the data set an image data set by 
implementing an optical tomography algorithm. The optical tomography algorithm may 
20 use factors related to determined probability distribution of photons attributable to the 
scattering character of the tissue being imaged. 

The controller may be arranged to activate the source and the detector to obtain a 
first selected distance between the input and detection ports, and the processor may be 
arranged to form the data set for the first distance. The processor may produce an image 
25 data set from the data set formed for the first distance. The controller may further be 
arranged to activate the source and the detector to obtain a second selected distance 
between the input and detection ports and is arranged to form another data set for the 
second distance. 



The processor may be arranged to correlate the first and second data sets by 



15 




The optical system may further include a display device constructed to receive 
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the image data set from the processor and to display an image. 

The optical system may further include a first oscillator and a phase detector. 
The first oscillator is constructed to generate a first carrier waveform at a first frequency 
on the order of 10 8 Hz, the first frequency having a time characteristic compatible with 
5 the time delay of photon migration from the input port to the detection port. The light 
source is coupled to the first oscillator and constructed to generate the light modulated 
by the first carrier waveform. The phase detector is constructed to determine change in 
waveform of the detected light relative to the waveform of the introduced light and 
measure therefrom the phase shift of the detected light at the wavelength, wherein the 
10 phase-shifted light is indicative of scattering or absorptive properties of the examined 
tissue region. The processor is arranged to form the data set based on the measured 
phase shift. This optical system may further include a second oscillator constructed to 
generate a second waveform at a second frequency. The detector is then arranged to 
receive a reference waveform at a reference frequency offset by a frequency on the order 
15 of 10 3 Hz from the first frequency and to produce a signal, at the offset frequency, 

corresponding to the detected radiation. The phase detector is adapted to compare, at the 
offset frequency, the detected radiation with the introduced radiation and to determine 
therefrom the phase shift. 

The optical system may further include an oscillator, a phase splitter, and first 
20 and second double balanced mixers. The oscillator is constructed to generate a first 
carrier waveform of a selected frequency compatible with time delay of photon 
migration from the input port to the detection port The light source is connected to 
receive from the oscillator the carrier waveform and is constructed to generate optical 
radiation modulated at the frequency. The phase splitter is connected to receive the 
25 carrier waveform from the oscillator and produce first and second reference phase 
signals of predefined substantially different phases. The first and second double 
balanced mixers are connected to receive from the phase splitter the first and second 
reference phase signals, respectively, and are connected to receive from the detector the 
detector signal and to produce therefrom a in-phase output signal and a quadrature 
30 output signal, respectively. The processor being connected to the double balanced 
mixers and arranged to receive the in-phase output signal and the quadrature output 
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signal and form therefrom the data set. 

The processor may be arranged to calculate a phase shift (@ x ) between the light 
introduced at the input port and the light detected at the detection port prior to forming 
the data set. 

The processor may arranged to calculate an average migration pathlength of 
photons scattered in the examined tissue between the optical input port and the optical 
detection port prior to forming the data set. 

The processor may further employ the pathlength in quantifying hemoglobin 
saturation (Y) of the examined tissue. 

The processor may be arranged to calculate a signal amplitude (A x ) determined 
as a square root of a sum of squares of the in-phase output signal and the quadrature 
output signal prior to forming the data set. 

The optical system may further include a narrow band detector connected to 
receive from the optical detector the detector signal and to produce a DC output signal 
therefrom. The processor then further determines a modulation index (M x ) as a ratio of 
values of the signal amplitude and the signal amplitude plus the DC output signal. 

The optical system may further include at least one oscillator connected to at 
least one light source. The oscillator is constructed to generate a carrier waveform of a 
selected frequency. The light source generate slight of a visible or infrared wavelength 
being intensity modulated at the frequency to achieve a known light pattern. The 
controller is constructed to control the emitted light intensity or phase relationship of 
patterns simultaneously introduced from multiple input ports, wherein the introduced 
patterns form resulting radiation that possesses a substantial gradient of photon density 
in at least one direction. This resulting radiation is scattered and absorbed over the 
migration paths. The detector is constructed and arranged to detect over time the 
resulting radiation that has migrated in the tissue to the detection port. The processor is 
further arranged to process signals of the detected resulting radiation in relation to the 
introduced radiation to create the data sets indicative of influence of the examined tissue 
upon the substantial gradient of photon density of the resulting radiation. 

The optical system may further include a phase detector constructed to detect the 
phase of the detected radiation and provide the phase to the processor. 
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The optical system may further include an amplitude detector constructed to 
detect the amplitude of the detected radiation and provide the amplitude to the processor. 

The phase relationship of light patterns introduced from two input ports may be 
180 degrees. 

The optical system may be constructed as described in U.S. Patent 5,1 19,815 or 
5, 386,827. This system includes a light source constructed to generate pulses of 
radiation of the wavelength, the pulses having a known pulse wave form of a duration on 
the order of a nanosecond or less. An optical detector is constructed to detect over time 
photons of modified pulses that have migrated in the tissue from the input ports. This 
system also includes an analyzer connected to the detector and adapted to determine a 
change in the pulse waveform shape of the detected pulses relative to the introduced 
pulses, at the employed wavelength. The processor then creates the data set based on the 
determined pulse waveform change. The processor may also be constructed and arranged 
to calculate the effective pafhlength of photons of the wavelength migrating between the 
input and detection ports in conjunction with creating the data set. The processor may 
also be constructed and arranged to calculate the scattering coefficient at the wavelength 
in conjunction with creating the image data set The processor may also be constructed 
and arranged to calculate the absorption coefficient at the wavelength in conjunction 
with creating the data set. 

The optical system may use the light source that produces relatively long light 
pulses and the processor that forms the data set by subtracting amplitude of two the 
pulses emitted from two input ports located symmetrically relative to one detection port. 

The optical system may be constructed to introduce and detect photons at two 
wavelengths selected to be sensitive to a tissue constituent. The tissue constituent may 
be an endogenous pigment or an exogenous pigment. The endogenous pigment may be 
hemoglobin. The exogenous pigment may be a selected contrast agent. 

Other advantages and features of the invention will be apparent from the 
following description of the preferred embodiment and from the claims. 



Brief Description of the Drawings 

Figs. 1 and 1 A show an optical module located on the right breast of a female 



subject. 

Figs. 2 and 2A show another embodiment of the optical module located on the 
right breast of the female subject. 

Figs. 3 and 3A show diagrammatically respective single wavelength and dual 
wavelength phase cancellation imaging systems that employ the optical module of Fig. 
1A or Fig. 2A. 

Fig. 3B is a timing diagram used by the imaging system of Figs. 3 and 3 A. 

Figs. 4 and 4A show diagrammatically another embodiment of the phase 
cancellation imaging system employing the optical module of Fig. 1 A or Fig. 2 A. 

Fig. 5 shows diagrammatically another embodiment of the phase cancellation 
imaging system employing the optical module of Fig. 1A or Fig. 2A. 

Fig. 6 shows schematically an amplitude cancellation imaging system using 
another embodiment of the optical module shown in Fig. 6A. 

Figs. 7, 7A and 7B show different embodiments of a cooling module used with a 
broad band light source, such as a tungsten light bulb. 

Fig. 8 shows diagrammatically another embodiment of the amplitude 
cancellation imaging system employing the optical module of Fig. 2A. 

Fig. 8A shows a circuit configuration for one element of the amplitude 
cancellation imaging system of Fig. 8. 

Fig. 8B is a timing diagram used by the imaging system of Fig. 8. 

Fig. 8C shows diagrammatically one channel of the amplitude cancellation 
imaging system of Fig. 8. 

Fig. 8D shows diagrammatically another embodiment of the amplitude 
cancellation imaging system of Fig. 8. 

Figs. 9A and 9B show respectively the left breast and the right breast during 
examination of a female patient lying supine face up. 

Figs. 10A, 10B, IOC, 10D, 11 A, 11B, 11C, 11D, 12A, 12B, 12C, 12D, 13A, 
13B, 13C, and 13D are back-projection images of breast tissue collected by the 
amplitude cancellation system of Fig. 8 

Fig. 14 is a two dimensional display of optical data from 30 studies of patients 
with suspicious masses. 



Figs. 15 and 15A are a two dimensional displays of optical data shown in Fig. 14 
using the contralateral breast data as a reference and using the breast model data as a 
reference, respectively. 

Figs. 16 A, 16B and 16C are histograms of data displayed in Figs. 14, 15 and 
5 15 A, respectively. 

Fig. 17 is a "four" dimensional graph summarizing 30 studies of patients with 
suspicious masses. 

Description of the Preferred Embodiments 

Referring to Figs. 1, 1A, 2 and 2A the right breast and the left breast of a female 

10 subject 8 are examined using an imaging system connected to an optical module 12 or 
14. Optical modules 12 and 14 include a multiplicity of light sources (e.g., laser diodes, 
LEDs, flashlight bulbs) providing light in the visible to infrared range and light detectors 
(e.g., photomultiplier tubes, Si diode detector, PIN, avalanche or other diode detectors), 
which may also include interference filters. The light sources and the light detectors are 

1 5 arranged to form selected geometrical patterns that provide a multiplicity of source- 
detector paths of photon migration inside the breast tissue. An optical examination 
system provides an in vivo optical data of the examined tissue, and the data may be 
processed to create an image. The image can show a location and size of an abnormal 
structure in the tissue, such as a tumor or bleeding. Furthermore, the optical data can 

20 provide a qualitative and quantitative measure (e.g., metabolic biochemistry, 

pathophysiology) of an abnormal tissue structure. (Alternatively, an optical module 
includes a multiplicity of optical fibers connected to one or several light sources and a 
multiplicity of optical detection fibers connected to one or several light detectors as 
described in the PCT Applications PCT/US96/00235 and PCT/US96/1 1 630 (filed 

25 January 2, 1 996 and July 12,1 996).) 

In one embodiment, optical module 12 includes nine laser diodes S l5 S 2 ,..., S 9 
and four photomultiplier tubes (PMTs) D,, D 2 , D 3 , D 4 . The laser diodes and PMTs are 
embedded in a pliable rubber-like material positioned in contact with the examined 
breast. There may be a Saran® wrap or similar material located between the laser diodes 

30 and the skin, and between the PMTs and the skin. Similarly, optical module 14 includes 
four laser diodes S,, S 2 , S 3 , S 4 and 27 silicon diode detectors D l5 D 2 , D 27 embedded in 



a pliable rubber-like material. The optical systems shown in Figs. 3 through 7 may be 
interfaced with optical module 12 or 14 for imaging of the breast tissue. Optical 
modules 12 and 14 have pairs of optical input ports symmetrically located (or 
equidistantly located) relative to an optical detection port, or have pairs of optical 

5 detection ports symmetrically located relative to an optical input port. In general, 

however, the ports do not have to be positioned symmetrically. The optical systems can 
vary the source or detector gain to account for any positional asymmetry or can 
introduce a selected asymmetry by adjusting the source or detector gain. 

Furthermore, the optical systems shown in Figs. 3 through 7, may be interfaced 

10 with two identical optical modules (12 or 14) located on the right breast and the left 
breast for lateralization, that is, comparative examination of the symmetric parts of the 
right breast and the left breast. For calibration, the optical module may also be placed on 
one or several breast models having the same scattering coefficient and the same 
absorption coefficient as the normal breast tissue of a female, for example, below 40 

15 years and above 40 years, and having different sizes and geometries. Alternatively, the 
breast tissue models may have other selected values of the scattering and absorption 
coefficient. 

Referring to Figs. 1A and 3, a phased array imaging system 15 is connected to 
optical module 12 with nine laser diodes S,, S 2 ,..., S 9 and four PMTs D l5 D 2 , D 3 , D 4 (e.g., 

20 Hamamatsu R928, Hamamatsu R1645u, T08) powered by a high voltage supply (not 
shown). Four laser diodes surround each PMT forming an equidistant arrangement (for 
example, different optical modules may use distances of 3.5, 7 and 10.5 cm). Switch 18 
connects laser diodes S„ S 2 ,..., S 9 to a phase splitter 20, which provides to the diodes an 
RF modulation signal having both a 0 degree phase and a 180 degree phase. Imaging 

25 system 1 5 also includes a 50 MHz single side band transmitter 22 connected by a phase 
lock loop 24 to a 50 MHz single side band receiver 26. Single side band (SSB) 
transmitter 22 is connected to a 1kHz oscillator 28, which provides a reference signal 30 
to a phase detector 32. SSB receiver 26 is connected to a switch 27, which connects one 
of the four PMTs (0.5 pV sensitivity) depending on control signals from a controller 19. 

30 The SSB transmitter-receiver pair can operate in the frequency region of 1 0-1 000 MHz 
(preferably 50-450 MHz). The phase noise of this apparatus is less than about 0.1°. The 



SSB receiver detects signal levels on the order of microvolts in a 2 KHz bandwidth. 
However, this narrow bandwidth limits the spread of switching of various light sources 
to approximately 1.0 msec, and thus the sequencing time for an entire image of 16 
source detector combinations can be ~ 1 sec. The system uses a 1 sec averaging time. 

5 Controller 19, connected to a personal computer (not shown), sequences laser 

diodes S,, S 2 ,..., S 9 so that two diodes receive 0° phase and 180° phase signals from 
splitter 20, every 0.1 sec. At the same time, controller 19 connects a symmetrically 
located PMT to SSB receiver 26. As shown in a timing diagram 40 (Fig. 3B), phased 
array imaging system 15 triggers two sources so that they emits modulated light of a 0° 

10 phase and a 1 80° phase for about 1 00 msec, and at the same time triggers a 

symmetrically located PMT. For example, when laser diodes 1 (S,) and 2 (S 2 ) emit light 
of the 0° and 180° phase respectively, detector 1 (D,) detects light that has migrated in 
the breast tissue. SSB receiver 26, which is phase locked with SSB transmitter 22, 
receives signal from detector 1 and provides output signal 34 to phase detector 32. 

15 Phase detector 32 measures the phase (36) of the detected light, and SSB receiver 26 
provides the amplitude (38) of the detected light. This phase detection circuit was— 
described in U.S. Patent 4,972,331, which is incorporated by reference. 

In the next cycle, controller 19 directs switch 18 to connect laser diodes 2 (S 2 ) 
and 3 (S 3 ), which emit modulated light of a 0° phase and a 180° phase, respectively, and 

20 detector 2 (D 2 ) detects light that has migrated in the breast tissue. Controller 1 9 also 
directs switch 27 to connect detector 2 to SSB receiver 26, which receives detection 
signal corresponding to the photons that have migrated from laser diodes 2 and 3 to 
detector 2. Again, phase detector 32 measures the phase (36) of the detected light, and 
SSB receiver 26 provides the amplitude (38) of the detected light. The duration of each 

25 pair of light flashes is 1 00 msec. The complete set of data for all source detector 

combinations is collected every 30 sec. A computer (not shown) stores the phase values 
and the amplitude values measured for the different combinations shown in timing 
diagram 40 and employs these values to create images of the examined tissue, as is 
described below. The computer uses the ADA2210 board for data acquisition. 

30 Before or after the above-described measurement, phased array imaging system 

1 5 may be calibrated on one or several models of the female breast. In the calibration 



procedure, the optical module is placed on the model and the imaging system collects 
the phase data and the amplitude data using the sequences shown in the timing diagram 
40. The imaging system may use different breast models having the same scattering 
coefficient and the same absorption coefficient as the normal breast tissue of a female, 

5 for example, below 40 years and above 40 years. Due to the nature of the visible or 
infrared optical radiation, the described optical imaging systems may be used for breast 
tissue examination of a woman of any age. Furthermore, the models may have different 
sizes and have a shape of a female breast during examination. For example, the model 
has a shape of the breast (or the breast and a portion of the chest) of a woman lying on 

1 0 her back during examination. 

Phased array imaging system 15 generates a "model" image for each wavelength 
employed. The model image may later be subtracted from the breast images to calibrate 
the system and also account for the boundary conditions of the light migrating in the 
tissue. Alternatively, phased array imaging system 15 is calibrated prior to taking 

15 measurement data and the gain on the light sources or the detectors is adjusted to obtain 
selected values. 

Referring to Figs. 1 A and 3 A, a dual wavelength phased array imaging system 45 
is connected to optical module 12 with nine 780 nm laser diodes S,, S 2 ,..., S 9 , nine 830 
nm laser diodes S la , S 2a ,..., S^, and the four PMTs D,, D 2 , D j5 and D 4 powered by a high 

20 voltage supply (not shown). Pairs of laser diodes S, and S u S 2 and S 2a , S 9 and S 9a 
are located next to each other and arranged to introduce modulated light at almost the 
same tissue locations. A switch 48 connects laser diodes S ,, S 2 ,..., S 9 to a phase splitter 
50, which provides to the laser diodes an RF modulation signal having both a 0 degree 
phase and a 180 degree phase. Similarly, a switch 48a connects laser diodes S la , S 2a ,..., 

25 to a phase splitter 50a, which provides to the laser diodes an RF modulation signal 
having both a 0 degree phase and a 1 80 degree phase. A 52 MHz SSB transmitter 52 is 
connected by a phase lock loop 54 to a 52 MHz SSB receiver 56, and a 50 MHz SSB 
transmitter 52a is connected by a phase lock loop 54a to a 50 MHz SSB receiver 56a. 
Both SSB transmitters 52 and 52a are connected to a 1kHz oscillator 58, which provides 

30 a reference signal 60 to phase detectors 62 and 62a. SSB receivers 56 and 56a are 
connected to one of the four PMTs by a switch 57 depending on control signals from 



controller 49. Controller 49, connected to a personal computer, sequences the laser 
diodes so that two pairs of the laser diodes receive 0° phase and 180° phase signals from 
splitters 50 and 50a, and at the same time controller 49 connects a symmetrically located 
detector to SSB receivers 56 and 56a. 
5 As shown in timing diagram 40 (Fig. 3B), phased array imaging system 45 

triggers for each wavelength two sources that emit modulated light of a 0° phase and a 
1 80° phase for about 100 msec and, at the same time, controller 49 connects the 
symmetrically located PMT. For example, switch 48 connects SSB transmitter 52 to 
780 nm laser diode 4 (S 4 ) to emit 52 MHz modulated light of a 180° phase and connects 
10 780 nm laser diode 5 (S 3 ) to emit 52 MHz modulated light of a 0° phase. At the same 
time, switch 48a connects SSB transmitter 52a to 830 nm laser diode 4a (S 4a ) to emit 50 
MHz modulated light of a 180° phase and connects 830 nm laser diode 5a (S 5a ) to emit 
52 MHz modulated light of a 0° phase. Simultaneously, switch 57 connects detector 1 
(Dj) to SSB receivers 56 and 56a to receive the detection signal corresponding to 
1 5 photons of both wavelengths that have migrated in the breast tissue. Phase detector 62 
provides the phase (66) of the detected 780 nm light, and phase detector 62a provides the 
phase (66a) of the detected 830 nm light for the selected geometry. Similarly, SSB 
receiver 56 measures the amplitude (68) of the detected 780 nm light and SSB receiver 
56a measures the amplitude (68a) of the detected 830 nm light. This operation is 
20 repeated for all combinations of sources and detectors shown in timing diagram 40. A 
computer (not shown) stores the phase values and the amplitude values (at each 
wavelength) measured for the different combinations shown in timing diagram 40. The 
computer then uses the measured values to create images using algorithms included the 
enclosed source code. 

25 Initially the system takes quick pictures to find the area of interest so that the 

optical module can be moved around to find an optimal geometry. Once found, the 780 
nm and 830 nm data (i.e., both the phase and amplitude data) is acquired and saved on a 
disk. 

Several phased array systems were described in the PCT application PCT/US 
30 93/05868 (published as WO 93/25 14 on December 23, 1 993), which is incorporated by 
reference. This PCT publication also describes the basic principles of phase and 



amplitude cancellation. The phased array imaging system uses a detector for detecting 
light emitted from equidistant sources located symmetrically with respect to the detector 
(or one source and several equidistant detectors located symmetrically). If two sources 
S, and S 2 emit modulated light having equal amplitude and a 0° phase and a 1 80° phase, 

5 detector D, located in the middle detects a null in the amplitude signal and detects a 
crossover between the 0° and 180° phase, i.e., a 90° phase, for substantially 
homogeneous tissue. That is, the detector is located on the null plane. In heterogeneous 
tissue, the null plane is displaced from the geometric midline. Nevertheless, the null 
establishes an extremely sensitive measure to perturbation by an absorber or scatterer. 

10 Furthermore, at the null condition, the system is relatively insensitive to amplitude 

fluctuations common to both light sources, and insensitive to inhomogeneities that affect 
a large tissue. The system has a high sensitivity to scattering provided that the scattering 
contrast is the same as the absorbing contrast. The system can readily observe shifts of 
50 to 60° of phase under altered blood volume or blood oxygenation conditions, where 

1 5 the phase noise is less than a 0.1° (s/n >400) for a 1 Hz bandwidth. The amplitude signal 
is little less useful in imaging since the position indication is somewhat ambiguous, i.e., 
an increase of signal is observed regardless of the displacement of the absorbing object 
with respect to the null plane, although this is remedied by further encoding of the 
sources. 

20 As described in the PCT application PCT/US 93/05868, the light sources excite a 

photon diffusion wave, which due to cancellation effects has a relatively long 
wavelength (-10 cm), determined by the scattering (u. s ' = 10 cm" 1 ) and absorption (u a = 
0.04 cm"') properties of the tissue. The photon diffusion wavelength of about 10 cm 
provides imaging in the "near field." The imaging system may use light sources of one 

25 or several optical wavelengths in the visible to infrared range, depending on the 

characteristic to be imaged (i.e., blood volume, blood oxygenation, a distribution of a 
contrast agent in the tissue, an absorbing constituent of the tissue, a fluorescing 
constituent of the tissue, or other). The phase signal at zero crossing detection is 
essentially a square wave "overloaded" signal. It is moderately insensitive to the 

30 changes of signal amplitude that may occur in imaging from proximal to distal source- 
detector pairs and is also moderately insensitive to ambient light. 



Referring to Fig. 4, in another embodiment, a phased array imaging system 100 
is used instead of imaging systems 15 or 45. Imaging system 100, connected to optical 
module 12 (shown in Fig. 1 A) having nine laser diodes S„ S 2 ,.„, S 9 and four PMTs D,, 
D 2 , D 3 , and D 4 , employs homodyne phase detection. A switch 1 02 connects laser diodes 

5 S,, S 2 ,..., S 9 to a phase splitter 1 04, which provides to the diodes an RF modulation 
signal having both a 0 degree phase and a 180 degree phase. Imaging system 100 also 
includes a 200 MHz oscillator 1 06 providing RF signal to a driver 108, which is 
connected to phase splitter 1 04. (Alternatively, an oscillator in the range of 1 0 - 1 000 
MHz, preferably 50 - 500 MHz, may be used.) A phase shifter 1 14 receives the drive 

10 signal (1 12) from driver 108 and provides the signal of a selected phase (e.g., a 0° phase 
change) to a 90° phase splitter 1 16. Phase splitter 116 provides a 0° phase signal (118) 
and a 90° phase signal (120) to double balance mixers (DBM) 122 and 124, respectively. 

A controller 140, connected to a personal computer, sequences laser diodes S u 
S 2 ,..., S 9 using switch 102 so that two diodes receive modulate signal at a 0° phase and a 

15 1 80° phase signals from splitter 104. At the same time, a controller 140 connects a 
symmetrically located PMT using a switch 130 to an amplifier 134. Amplifier 134- 
provides a detection signal (136) to double balance mixers 122 and 124, and to a DC 
detector 138. Double balance mixer 122 receives the detection signal (136) and the 0° 
phase reference signal (118) and provides an in-phase signal I (144). Double balance 

20 mixer 124 receives the detection signal (136) and the 90° phase reference signal (120) 
and provides a quadrature signal R (142). DC detector 138 provides DC signal (146). 
The in-phase signal I and quadrature signal R specify the phase (8 = tan"'I/R) of the 
detected optical radiation and the amplitude (A = (R 2 + P)' 1/2 ) of the detected optical 
radiation. This phase detection circuit was described in U.S. Patent 5,553,614, which is 

25 incorporated by reference. 

Similarly as for imaging systems 15 and 45, imaging system 100 directs 
controller 140 to sequence the laser diodes and the PMT detectors using timing diagram 
40. The computer stores the phase value and the amplitude value measured for each of 
the combinations and generates images described below. 

30 Fig. 4A shows diagrammatically one portion of phase cancellation, phased array 

imaging system 100. The depicted portion of imaging system 100 includes two laser 



diodes LD,, and LD 2 and a light detector D,, which are included in optical module 12 or 
14. Oscillator 106 provides carrier waveform having a frequency in range of 30 to 140 
MHz. The carrier waveform frequency is selected depending on the operation of the 
system. When time multiplexing the light sources using switch 102, then the carrier 
waveform is modulated at a lower frequency, e.g., 30 MHz to afford switching time. 

When no time multiplexing is performed, oscillator 1 06 operates in the 1 00 MHz 
region. Splitter 1 04 splits the oscillator waveform into 0° and 1 80° signals that are then 
attenuated by digitally controlled attenuators 107 A and 107B by 0% to 10% in 
amplitude. The phase of the attenuated signals is appropriately shifted by digitally 
controlled phase shifters 109A and 109B in the range of 10° - 30°' and preferably 20° in 
phase. Laser drivers 108 A and 108B drive LD , and LD 2 , respectively, which emit light 
of the same wavelength, for example, 780 or 800 nm. After the introduced light 
migrates in the examined tissued, a PMT detector D, amplifies the detected signals 
having initially the 0 and 180° phases. As described above, for homogeneous tissue and 
symmetric locations of LD ]; LD 2 and D„ the output of the PMT is 90°, i.e., halfway 
between 0° and 180° and the amplitude is close to zero. The personal computer (PC) 
adjusts the attenuation provided by attenuator 107B and the phase shift provided by 
phase snifter 109B so that detector D, detects phase nominally around 25° and 
amplitude nominally around £ 10 millivolts for homogeneous tissue. This signal is 
connected to amplifier 134 and to the IQ circuit 139. The cosine and sine signals are fed 
into the personal computer, which takes the amplitude (the square root of the sum of the 
squares of I and Q) and the phase angle (the angle whose tangent is I/Q) to give outputs 
of phase around 25° and amplitude signals around 10 millivolts. The personal computer 
also adjusts the reference signal to the IQ to have the phase (j> 3 between 10° to 30° and 
preferably around 25°, i.e., phase shifter 114 provides to the IQ circuit 139 the reference 
phase having a value selected by the combination of phase shifters 109A and 109B. 

In a currently preferred embodiment, splitter 104 is a two way 180° power 
splitter model number ZSCJ-2 1, available from Mini-Circuits (P.O. Box 350186, 
Brooklyn, New York 1 1235-0003). The phase shifters 109 A, 109B and 1 14 and 
attenuators 107 A, and 107B are also available from Mini-Circuits, wherein the 
attenuators can be high isolation amplifier MAN- IAD. IQ demodulator 139 is a 



demodulator MIQY-140D also available from Mini-Circuits. 

The system obtains the initial values of attenuator 107B (A 2 ) and phase shifter 
109B (4> 2 ) on a model or a symmetric tissue region (e.g., the contralateral breast that is 
tumor free). The entire probe is calibrated on a tissue model by storing the calibration 
values of A 2 and (j> 2 for the various source-detector combinations (i.e., the baseline 
image). The probe is then moved to the breast, for example, and the phases and 
amplitudes are detected for the various source and detector combinations. When the 
contralateral tumor free breast is used as a model, the probe is transferred to the 
contralateral breast (taking note to rotate the probe because of the mirror image nature of 
the breast physiology) and then the images are read out from all the source-detector 
combinations to acquire the tissue image. There is no limitation on multiplexing as long 
as the bandwidth of Ft andF 2 is recognized as being the limiting condition in the 
system normalization. It should be noted that normalization must be accurate and 
without "dither" and therefore, a significant amount of filtering in F, and F 2 , i.e., less 
than 10 Hz bandwidth. If <|> 2 is adjusted over a large range, there will be an amplitude- 
phase crosstalk. Thus, the system may adjust phase and then amplitude and repeat these 
adjustments iteratively because of the amplitude phase crosstalk. The control of A , and 
cb, provides even a greater range of control, where obviously inverse signals would be 
applied to them, i.e., as the Ai(b, signals are increased, the A 2 , <j> 2 signals would be 
decreased. Both A 2 and 4> 2 can be controlled by PIN diodes, to achieve an extremely 
wideband frequency range. However, since signal processing controls the bandwidth of 
the feedback system, that either PIN diode or relay control of the phase and amplitude is 
feasible for automatic compensation. If, in addition, dual wavelength or triple 
wavelength sources are used, each one of them must be separately calibrated because no 
two light sources can be in the same position relative to the imaged tissue (unless, of 
course, they are combined with optical fibers). 

Referring to Fig. 5, in another embodiment, a dual wavelength phased array 
imaging system 150 is used instead of imaging systems 15, 45 or 100. Imaging system 
150, connected to optical module 12 (shown in Fig. 1 A) having nine 760 nm laser diodes 
S„ S 2 ,..., S 9 nine 840 nm laser diodes S !a , S 2a ,..., S 9a and four PMTs D„ D 2 , D 3 , and D 4 is 
based on heterodyne phase detection. A switch 152 connects the laser diodes to a phase 



splitter 1 54, which provides to the diodes an RF modulation signal having both a 0 
degree phase and a 1 80 degree phase. Imaging system 150 employs a mixer 1 65 
connected to a 200 MHz oscillator 160 and 200.025 MHz oscillator 162 (Alternatively, 
oscillators operating in the range of 10 - 1000 MHz, preferably 50 - 500 MHz, may be 
used.) Mixer 165 provides a 25 kHz reference signal (168) to an adjustable gain 
controller 177. Oscillator 162 connected to power amplifier 163 provides a 200.025 
MHz reference signal (170) to the second dynode of each PMT detector for heterodyne 
detection. Each PMT detector provides a 25 kHz detection signal (172) to a switch 178, 
which in turn provides the signal to a 25 kHz filter 1 80. A phase detector 1 84 is 
connected to an adjustable gain controller 182, which provides a filtered and amplified 
detection signal (186) and to adjustable gain controller 177, which provides the reference 
signal (188). Phase detector 184, connected to a switch 190, provides the detected phase 
value for each wavelength. This phase detection circuit was described in U.S. Patent 
5,1 87,672, which is incorporated by reference. Another type of phase detection circuit 
was described in U.S. Patent 5,564,417, which is incorporated by reference. 

Similarly as described above, controller 175, connected to a personal computer, 
sequences laser diodes S„ S 2 ,..., S 9 or laser diodes S la , S 2a ,..., S 9a using switch 152 so that 
two diodes emitting the same wavelength receive 0° phase and 180° phase signals from 
splitter 154. At the same time, controller 175 connects a symmetrically located PMT 
using a switch 178 to filter 180 and adjustable gain controller 182. Phase detector 184 
provides the measured phase. Imaging system employs timing diagram 40 (Fig. 3B); 
however, since the two wavelength light is not frequency encoded, laser diodes S l5 S 2 ,..., 
So or laser diodes S la , S 2a ,..., S 9a are triggered in each sequence. The computer stores the 
phase values measured for the different combinations and generates images described 
below. 

Referring to Fig. 6, in another embodiment, an amplitude cancellation imaging 
system 200 uses an optical module 212 shown in Fig. 6B. Optical module 212 includes 
twelve light sources SI, S2, SI 2 and four light detectors Dl, D2, D3, and D4 
mounted on a plastic or rubber foam material. The light sources and the light detectors 
) are located on a geometrical pattern that provides sixteen source-detector combinations 
(CI, C2, C16) having a selected source-detector separation. The separation may be 



2.5 cm to produce about 1 .25 cm average light penetration. (Several modules with 
different source-detector separations may be used to obtain several two dimensional 
images of different tissue depths. Alternatively, a single module may include source 
detector combinations providing different separations.) The light sources are 1 W 
tungsten light bulbs, which emit broad band non-modulated light. The light detectors 
are silicon diodes, each equipped with an interference filter transmitting a 10 ran wide 
band centered at 760 nm and 850 nra. The 760 nm and 850 nm wavelengths are selected 
to detect oxyhemoglobin and deoxyhemoglobin in the examined tissue. 

Optical module 212 is connected to an analog circuit 202, which includes a 
source circuit 204 for controlling sources SI, S2, ... S12. Optical module 212 is 
connected to a detector circuit 206, which controls diode detectors Dl, D2, D3 and D4. 
In general, imaging system 200 can turn ON each source for a selected period in the 
range of 10" 6 sec. to 0.1 sec, and one or several symmetrically located detectors are 
turned on simultaneously or sequentially to collect optical data. Specifically, as 
provided in Appendix B, one of sources SI, S2, ... S12 is turned ON for 500 msec and 
the emitted light is introduced into the tissue from the corresponding input port. The 
introduced photons migrate over banana shaped paths in the examined tissue to a 
detection port. The corresponding detector is triggered 200 msec, after the source and 
collects light for 200 msec. A detection circuit 206 receives a detector signal from the 
diode detector. Detection circuit 206 enables correction for the dark current/noise that 
comprises background light, DC offset of the operational amplifiers, photodiode dark 
current, temperature effects on the outputs of individual components and variations due 
to changing environment. 

Imaging system 200 performs data acquisition in four steps synchronized by its 
internal oscillator. The first step is performed by having the light sources OFF. The 
detector output is directed to an integrator 216 and integration capacitor 218 is charged 
to the dark level voltage. In the second step, the light source is turned ON and after 200 
msec the preamplifier output that corresponds to the intensity of the detected light is 
directed to integrator 216 in a way to charge capacitor 21 8 with current of polarity 
i opposite to the polarity of the charging current in the first step. This is achieved using 
an appropriate ON/OFF combination of switches A and B. The voltage of capacitor 218 



is charging to a value that, after 200 msec, represents the total detected intensity minus 
the dark level noise signal. In the third step, both switches A and B are turned OFF to 
disconnect both the positive unity gain and the negative unity gain operational amplifiers 
(220 and 222). Then, the output of integrator 21 8 is moved via switch C to an analog-to- 
5 digital converter and the digital signal is stored in the memory of a computer. In the 
fourth step, the switches A, B and C are open and switch D is closed in order to 
discharge capacitor 21 8 through a 47K resistor. At this point, the circuit of integrator 
216 is reset to zero and ready for the first step of the detection cycle. 

Alternatively, analog circuit 202 may be replaced by a computer with an analog- 
1 o to-digital converter and appropriate software that controls the entire operation of optical 
module 212. An algorithm that controls the sources and the detectors of optical module 
212 in a similar way as described above. The detected dark level noise signal is digitally 
subtracted from the detected intensity of the introduced light. 

The collected data sets are processed using an imaging algorithm. The imaging 
1 5 algorithm calculates the blood volume of the examined tissue for each source-detector 
combination for each data set. The imaging algorithm can also calculate the 
oxygenation of the examined tissue for each source-detector combination. 

The blood volume or oxygenation images can be subtracted from "model" 
images. The blood volume image can be subtracted from the oxygenation image to 
20 create congruence data to localize and characterize a tissue anomaly. That is, the 
imaging algorithm creates an image using the differential image data sets. Prior to 
creating the image, an interpolation algorithm is employed to expand the differential 
image data set, containing 16 (4x4) data points, to an imaging data set containing 32x32 
image points. 

25 Alternatively, the computer uses a back-projection algorithm known in computed 

tomography (CT), which is modified for light diffusion and refraction and the banana 
like geometry employed by the optical imaging system. In the optical back-projection 
algorithm, the probabilistic concept of the "photon migration density" replaces the linear 
relationship of ballistically transmitted X-rays, for the beam representing pixels. The 

30 photon migration density denotes a probability that a photon introduced at the input port 
will occupy a specific pixel and reach the detection port. For different types of tissue, 



the phase modulation spectrophotometer provides the values of the scattering and 
absorption coefficients employed in the probability calculations. In the image 
reconstruction program, the probability is translated into a weight factor, when it is used 
to process back-projection. Back-projection algorithms known in CT may be used. The 
5 back-projection averages out the values of information that each beam carries with the 
weighting in each pixel. A weighting algorithm for creating a photon density image may 
be used in the back-projection reconstruction algorithm described above. 

A method for correcting blurring and refraction used in the back-projection 
algorithm was described by S.B. Colak, H.Schomberg, G.W.'t Hooft, M.B. van der Mark 
1 o on March 12, 1 996, in "Optical Back-projection Tomography in Heterogeneous 

Diffusive Media" which is incorporated by reference as if fully set forth herein. The 
references cited in this publication provide further information about the optical back- 
projection tomography and are incorporated by reference as if fully set forth herein. 
Another embodiment of the amplitude cancellation imaging system 200 uses 
15 optical module 14 shown in Fig. 2A. In this arrangement, four centrally located light 
sources SI, S2, S3, and S4 and 21 detectors Dl, D2, D21 provide a multiplicity of 
symmetric photon migration paths for each source. For example, source SI is turned 
ON for a period in the range of 10" 6 sec. to 0.1 sec. The source emits non-modulated 
light into the breast tissue. Symmetrically located detectors Dl and Dl 1 are ON 
20 simultaneously to collect introduced photons migrating over substantially symmetric 
paths. For normal breast tissue, detectors Dl and Dl 1 detect light of the same intensity, 
and thus the differential signal is zero, i.e., the detected amplitude are canceled. Imaging 
system 200 collects the differential data for a multiplicity of symmetric photon migration 
paths and generates an image of the examined tissue. Imaging system 200 may collect 
25 optical data for several wavelengths and generate blood volume images and blood 

oxygenation images for the examined tissue. Amplitude cancellation imaging system 
200 may also use a second identical optical module 14 placed on the opposite breast. 
The blood volume or oxygenation images collected for the two breasts may be 
subtracted to provide a differential image, which will further emphasize a tissue 
30 abnormality located in one breast. 

Alternatively, the amplitude cancellation imaging system uses light modulated at 



frequencies in the range of 0.1 to 100 kHz. The system employs the above-described 
algorithm, but the light sources emit frequency modulated light and the detectors, each 
connected to a lock-in amplifier, detect light modulated at the same frequency. This 
lock-in detection may further increase the signal to noise ratio by eliminating external 
5 noise. The detected light intensities are processed the same way as described above to 
image the examined tissue. 

Figs. 7, 7A and 7B show different embodiments of a cooling module used with a 
broad band light source or light guides, where they are positioned close to the skin. In 
this arrangement, there is trapped heat that frequently causes an uncomfortable 
10 temperature. Fig. 7 depicts a cooling module 230, which surrounds light sources 232 A 
and 232B. Cooling module 230 includes a fan 234 and a set of air passages 236. In a 
similar design, two fans are juxtaposed on each side of one or more light bulbs to form 
an "open frame" so that the fans blow not only upon the light sources, but upon the skin 
itself. The cooling module enables a power increase on the light sources, but no increase 
1 5 of heat upon the skin itself, which remains under comfortable conditions. 

Fig. 7A depicts a cooling module 240 for cooling light guides. Light guides 242 
deliver light and heat to the skin. A cooling ring 244 includes an air inlet 246 and a set 
of air passages 248 (or jets) for providing air flow to the irradiation location. Fig. 7B 
depicts a cooling module 250 constructed to air cool a light barrier 252. Light barrier 
20 252 has similar optical properties as the light barrier described in the PCT application 
PCT/US92/04153 (published on November 26, 1992 as WO 92/20273), which is 
incorporated by reference. This embodiment utilizes the advantages of the light barrier 
and enables the use of higher light intensities. Cooling module 250 includes air inlets 
252A and 252B, which provide air to a set of conduits and openings that deliver air to 
25 the skin near light source 254. Compressed air may also be used. 

The intensity regulations for delivering continuous otherwise noncoherent light 
to the skin often depend on the temperature rise of the skin itself. For examination of 
large tissue volumes or deep tissues (i.e., where there is a large separation between the 
optical input and optical detection ports) relatively large light intensities are needed. 
30 Under conditions of prolonged even low level illumination, the skin may become 

uncomfortably warm and may blister. However, the erythemic effects are much smaller 



in the NIR, where the delivered heat is a factor, than they are in UVA and UVB, where 
cancer-producing damage may occur (but is not known for the NIR). The effect of the 
cooling air is not just convection of warm air away from the skin, but it enhances the 
evaporation of perspiration from the skin. Thus, as soon as the skin temperature rises 
5 and perspiration is initiated, greatly enhanced cooling is obtained with the forced air 
increasing the evaporation. 

Referring to Fig. 8, an amplitude cancellation imaging system 260 is used instead 
of imaging systems 15, 45, 100, 150, or 202. Dual wavelength amplitude cancellation 
imaging system 260 is connected to optical module 14 shown in Fig. 2 A and includes 
10 four 750 nm laser diodes S,, S 2 , S 3 , and S 4 , four 830 nm laser diodes S Ia , S 2a , S 3a , and 
S 4a , and twenty-one silicon diode detectors D,, D 2 , D 21 . Each detector is connected to 
a preamplifier and an adjustable gain controller that may be used initially for calibration. 
The detector outputs are switched by a switch 262 by a controller 264 so that analog-to- 
digital converters 266 and 266a receive 750 nm and 830 nm data, respectively, from two 
15 symmetrically located detectors. A computer 270 stores the detected values measured 
for the different combinations using algorithms employed in the enclosed source code. 
The computer also generates images described below. Another type of amplitude 
detection circuit was described in Figures 1 1 through 13 and the corresponding 
specification of U.S. Patent 5,673,701, which is incorporated by reference as if fully set 
20 forth herein. 

Also referring to Figs. 8A and 8B, the controller sequences an oscillator 261 so 
that each source emits a 50 jisec light pulse as shown in timing diagram 272. The 
system sequences through the various source/detector combinations in approximately 
one msec, and averages the imaged data over 8 sec to get a very high signal to noise 

25 ratio. The circuit configuration for one element of imaging system 260, i.e., 754 nm 
sources S,, S 2 and 830 nm sources S la , S 2a , and two symmetrically positioned detectors 
D 3 and D n , is shown in Fig. 8A. The light intensities detected for the symmetrical 
locations are subtracted in a digital or analog way. The computer stores all data detected 
for the two wavelengths for generating tissue images. 

30 Fig. 8C shows diagramatically a single channel 260A of the time multiplex 

imaging system 260. Detector D t detects light emitted from light source S, emitting 



light pulses of the duration of about 50 usee. The detector signal is amplified and 
provided to a sample-and-hold circuit and filter. Detector D, is a silicon diode detector 
that has the detection area of about 4X4 mm and includes a pre-amplifier. The filtered 
signal 272 is provided to an AGC 274, which adjusts the amplitude of the signal based 

5 on a control signal from a personal computer. The personal computer has normalization 
amplitudes for the individual source-detector combinations. 

Amplitude cancellation imaging system 260 is normalized on a tissue model by 
detecting signals for the individual source-detector combinations and appropriately 
normalizing the detected signal using the AGC control. The individual 

10 normalization/calibration amplitudes form a baseline image that is stored in the 

computer. As described above, the baseline image may also be acquired on a symmetric 
tissue region, such as the contralateral breast for breast tissue examination or the 
contralateral brain hemisphere for brain tissue examination. The normalization process 
can be repeated several times to account for drifts in the individual elements. During the 

1 5 measurement process, the personal computer can adjust the gain of each AGC 3 14 based 
on the calibration values that account only for the electronic drift. Then, the defected 
image is subtracted from the baseline image of the examined tissue. Alternatively, while 
collecting the measurement data on the examined tissue, the measurement image is 
subtracted from the baseline image to create the tissue image that includes any tissue in 

20 homogeneities such as a tumor or bleeding. The sample-and-hold circuit maybe an 

analog circuit or the sample-and-hold function, including the filtering, may be performed 
digitally. 

Fig. 8D shows diagramatically an amplitude cancellation imaging system 
employing a frequency multiplex method. Amplitude cancellation system 300 includes 

25 21 oscillators 302 operating a frequencies in the range of 1 kHz to 100 kHz. Each 

oscillator 302 drives a light source 304 (for example, a laser diode or LED), which emits 
an intensity modulated light into the examined tissue. Each light detector 306 (for 
example, a photomultiplier, an avalanche photodiode PIN detector or a silicon detector) 
detects the intensity modulated light and provides a detector signal to an amplifier 308. 

30 The amplified detector signal is provided to a processing channel 310, which includes a 
band pass filter 312, an AGC 3 14, a lock-in amplifier 3 1 6, and a filter 318. Filter 312 



filters the detector signal, and AGC 314 adjusts the amplitude according to the input 
signal from a personal computer. Lock-in amplifier 3 1 6 receives the amplified signal 
315 and a reference signal 320 from oscillator 302. Lock-in amplifier 3 1 2 provides 
amplitude signal 3 1 7 to filter 318. Processing channel 310 may be an analog channel or 

5 a digital channel. 

In the amplitude cancellation system 310, all light sources emit light at the same 
time into a selected tissue region. Each light source is modulated at a distinct frequency 
in the range of 1 kHz to 1 00 kHz. In order to resolve the modulated light signals and 
attribute them to the individual light sources, the oscillators operate at frequencies 1 

10 kHz, 2 kHz, 4 kHz, 8 kHz, 16 kHz, ... Filters 312 and 3 1 8 are designed to provide only 
the detection signal from a selected light source, and lock-in amplifier 312 provides the 
amplitude of the signal at the selected frequency. Frequency multiplex system 300 is 
calibrated the same way as the time multiplex system 260, and the 
normalization/calibration amplitude values are also stored in the personal computer. The 

15 images are processed as described above. 

All above-described imagers will achieve a higher spacial resolution of the 
imaged tissue by increasing the number of sources and detectors. Furthermore, the 
sources and detectors may form various 1 dimensional, 1 .5 dimensional, or 2 
dimensional arrays as described in the above-referenced documents. 

20 In a clinical study, which is provided only as an illustrative example, a number of 

female patients were examined using three different positions. First, the patients were 
sitting upright supporting the breast by the imager optical module. Second, the patients 
were lying face down with the breast on the optical module pad, and third, the patients 
were lying supine face up with the breast spread over the chest as evenly as possible. 

25 Third, the patients were lying on their back. This position with the flattened adipose 
breast is the optimal position for imaging of a suspicious mass in patients over 40 years 
old. The third position is also preferred for patients below 40 years since the matching 
tissue background of the chest muscles and the ribcage provides suitable boundary 
conditions. In this position, no matching fluid was used. The thickness of the examined 

30 tissue with this geometry is usually no more than 4 to 5 cm even for a large breast. 
Initially, the imager was first calibrated on a breast model. During the 
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examination, the patient or the attendant held optical probe 14 over a designated portion 
of the breast. The mirror image region on the contralateral breast was also recorded. 
The images were acquired in 8 sec from each breast taking advantage of a priori 
information obtained by palpation, X-ray mammography, ultrasound examination, JMKL, 
5 or a needle biopsy. Other images were taken having only palpation information so that 
the 10 x 10 cm 2 active area of optical probe 14 covered the suspicious mass. The 
positioned probe may or may have not included the nipple, which is clearly visible on 
the optical images. 

The optical images were created using a back-projection algorithm. The images 

10 may be displayed in the format of the left breast data minus the model data, the right 
breast data minus the model data for each wavelength (e.g., 750 and 830 nm). 
Alternatively, the model calibration may be performed by adjusting the detector gains 
prior to the breast measurements. Furthermore, the images may be the differential 
between the right breast and the left breast for each wavelength to emphasize any tissue 

1 5 difference, such as a suspicious mass, which is unlikely located symmetrically in both 
breasts. The images may also be processed to image blood volume and blood 
oxygenation of the examined tissue of each breast. The blood volume image is the sum 
of 0.3 times the 750 nm data and 1.0 times the 830 nm data. The blood deoxygenation 
image is the difference of the 750 nm and the 830 nm data. The above coefficients were 

20 derived from blood tests in model systems. The images have the highest specificity and 
sensitivity when the contralateral breast data is used as a baseline and both the blood 
volume data and the hemoglobin deoxygenation data is imaged and positional compared. 

The blood volume and hemoglobin deoxygenation images provide an important 
tool in characterizing a suspicious anomaly in the examined breast. While the blood 

25 volume, hemoglobin oxygenation images and hemoglobin deoxygenation images, as 
well as the single wavelength images, are useful in locating an abnormal tissue region 
(i.e., detecting the abnormal structure), these images are also used to characterize the 
metabolism or pathology of the suspicious tissue anomaly. Specifically, an increased 
blood volume signal is observed with respect to the adipose tissue background due to the 

30 increased vascularity of a tumor as a consequence of angiogenetic factors. These factors 
include actively metabolizing regions and necrotic/apoptotic regions of the tumor. On 



the other hand, the hemoglobin deoxygenation signal is related to metabolic intensity. 
That is, the balance between oxygen delivery and oxygen uptake, which in tumors is 
usually balanced in favor of oxygen uptake exceeding oxygen delivery. The increased 
oxygen uptake occurs particularly for those tumors that are aggressively growing, and 
5 may well be metastatic. 

By selecting an appropriate wavelength, or several wavelengths, sensitive to an 
optically active tissue property, the imaging system can non-invasively characterize a 
tissue anomaly. The above-mentioned wavelengths are sensitive to hemoglobin and 
hemoglobin oxygenation, but other wavelengths sensitive to absorption by any tissue 
10 constituent may be used. Furthermore, an optical contrast agent (e.g., cardiogreen, 
indocyanine green) may be injected intravenously. The imaging system will then use a 
wavelength sensitive to the administered contrast agent. The regions of increased blood 
volume will also have a higher content of the contrast agent. 

Alternatively, differences in tissue scattering may be imaged. Due to differences 
15 in the optical refractive index, different types of tissue and different tissue solutes scatter 
light differently. The above-described imaging systems are also sensitive to scattering 
changes. The imaging system may use a wavelength that does not exhibit absorption 
changes for different types of tissue and different tissue solutes, but exhibits differences 
in scattering. 

20 The following are examples of clinical data acquired by the amplitude 

cancellation system of Fig. 8 employing optical module 14. Figs. 9A and 9B show 
respectively the left breast and the right breast during examination with a patient lying 
supine face up (i.e., the third position). The right breast (Fig. 9B) included a suspicious 
area 301 located below the nipple 303. Figs. 10A through 10D are back-projection 

25 images collected by the amplitude cancellation system of Fig. 8 employing optical 
module 14 as shown in Figs. 9A and 9B. Figs. 10A and 10B are images of the right 
breast corrected for the model data at 750 nm and 830 nm, respectively. The x and y 
coordinates are in millimeters and the intensity in 0 V to 3 V with respect to the model. 
Fig. 10C is a blood deoxygenation image of the right breast corrected for the model data. 

30 Fig. 1 0D is a blood volume image of the right breast corrected for the model data. 

Figs. 1 1 A and 1 IB are images of the right breast data minus the left breast data 
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at 750 nm and 830 nm, respectively (No. 18). Fig. 1 1C is a blood deoxygenation image 
of the right breast data minus the left breast data. Fig. 1 ID is a blood volume image of 
the right breast minus the left breast data. These images have a high signal level and 
exhibit high congruence of the tumor. 
5 Figs. 12A through 12D show optical images of a patient who had a "suspicious 

area" identified below the nipple of the left breast by X-ray mammography. Figs. 12A 
and 12B are images of the right breast data minus the left breast data at 750 nm and 830 
nm, respectively. Fig. 12C is a blood deoxygenation image of the right breast data 
minus the left breast data. Fig. 1 2D is a blood volume image of the right breast minus 
10 the left breast data. Figs. 12C and 12D exhibit only a blood volume increase of 0.1 8 V 
and a small negative value for deoxygenation of about 0.4. Therefore, these optical 
images would not be considered abnormal. 

Figs. 13 A and 13B are images of the right breast corrected for the model data at 
750 nm and 830 nm, respectively. Fig. 13C is a blood deoxygenation image of the right 
15 breast corrected for the model data. Fig. 1 3D is a blood volume image of the right breast 
corrected for the model data. 
Discussion of the Images 

Fig. 14 is a two dimensional plot of 30 studies of patients with suspicious masses 
based on mammographic data guidance provided as an example. Specifically, Fig. 14 
20 displays angiogenicity and hypermetabolic activity of the measured masses as the blood 
volume increases and the deoxygenation increases. All patients had a palpable mass, 
and/or an X-ray mammography and/or ultrasound detected suspicious breast mass. 
Patients nos.13, 16, 18 and 24 were identified as having a cancerous mass, while patients 
17, 20, 25 and 28 had a significant risk of having a cancerous mass. 
25 The two-dimensional nature of Fig. 14 permits the assignment of sensitivity and 

specificity according to the signal strength of blood volume data and the signal strength 
of deoxygenation data. I divided the region of signals in Fig. 14 into four zones. Zone I 
was defined for blood volume values above about 2.4 V and deoxygenation values above 
about 1 .4 V. Zone II, located below Zone I, was defined for blood volume values above 
30 about 1.7 V and deoxygenation values above about 0.75 V. Zone III, located below 
Zone II, was defined for blood volume values above about 1.3 V and deoxygenation 



values above about 0.2 V. Zone IV was located below Zone III. In Zone I three out of 
four patients had a cancerous mass. In Zones III and II three out of five patients had a 
cancerous mass (not counting patient 30). 
Bilateral Signals. 

5 In these cases, bilateral congruence affords a diagnostic aid. In some cases, 

particularly in cancer bearing breast (point 13 (GR)), the left and right breast signals 
measured with respect to the model showed similar shapes and thus subtraction of the 
left and right mammographically dense breasts failed to give a recognizable object in 
the designated outer, lower portion of the left breast. While indeed a residual shape was 

10 seen at the bottom of the image, the region was not intense; the blood volume and the 
deoxygenation signals were about 0.9 and 0.4 V. However, when the model is used as a 
reference, the middle, left portion of the image gives signals 2 and 0.9 V. In this case, 
the signal is placed in Zone II in Fig. 14. However, a further study of such masses is 
desirable. 

15 Subject No. 30 (likely normal) showed distinctive signals in the left portion of 

the right breast, but showed no distinctive signals in the region of the right breast, which 
was heterogeneous and dense. This subject (No. 30) showed no significant features in 
the data of the left breast minus model, but the data of the right breast minus model 
showed a clear anomaly in both blood volume and deoxygenation. This anomaly 

20 appeared on the data of the right breast on the inner side along the midline providing a 
highly significant signal from a congruent volumes. This subject was scored at 2.5 V and 
1 .9 V having heterogeneously dense breast (43 year old) but no mammographically 
detected tumor. 

Studies referenced to the model breast . 

25 No. 25 (KS) was found to have a mass near the nipple, 9-11 o'clock in the left 

breast. Both breasts were characterized by significant blood volume and deoxygenation 
signals as if the breast had unusual glandular or active metabolic characteristics. 
However, there was a very clear spot at the mammographically indicated region in the 
deoxygenation image with respect to the models, but little signal in the congruent 

30 position in the blood volume minus model image. On subtraction two, there is poor 

congruence, again a small image in deoxygenation, and large non-congruent image at the 
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mammographic indication. Thus, I scored the optical data of subject 25 against the 
model (i.e., the left breast minus model) because of background interference from the 
contralateral breast. Subject 15 lacks the optical data on the model, but a shape in the 
blood volume with no congruence in oxygenation was detected. 

5 The optical images of subject 8 (LA) had a "suspicious mass" below the nipple in 

both left breast and right breast. In fact, there was rather good congruence between the 
left and right breasts in the blood volume images but the "oxygenation mass" was distant 
from the "suspicious mass". Thus, in the difference blood volume image the mass 
disappeared with no image left in the suspicious region. 

10 A 23 year old subject (No. 9) showed a mammographically designated lower 

middle mass in the left breast. There is a congruent blood volume middle breast signal in 
both left and right breast showing a wide diffuse area of low intensity of deoxygenation 
in both breasts (about 1 .8 V). The subject had similar anomalies in both breasts, with 
congruent deoxygenation. This subject was scored on the left breast only providing the 

15 values of about 1 . 8 and 0.5V, respectively . 

A 27 year old patient (No. 28) with fibroadenoma in the left breast near the 
midline displayed a very large blood volume signal of 3.8 V and a relatively small and 
diffuse deoxygenation signal with very poor congruence (0.7 V) placing the signal in 
Zone III 

20 Nipple Signal 

In the contralateral (the right breast), the nipple shows as about a 1 cm shape in 
both the blood volume and deoxygenation images proving low signal levels of about 1 
and 1.2 V, respectively (see Figs. 13 A through 13D). 
Two Dimensional Charts 

25 Figs. 14 through 17 are provided to display and analyze the optical data. Fig. 14 

displays the optical data evaluated based on mammographic data. Fig. 15 shows a chart 
of optical data displayed in Fig. 14 based upon "suspicious mass" guidance only, using 
the contralateral breast data as a reference. Fig. 15A shows a chart of optical data 
displayed in Fig. 1 4 based, however, upon "suspicious mass" guidance only, using the 

30 breast model data as a reference. The difference between Fig. 14 and Figs. 15 and 15A 
is that in Fig. 14 the calculated data is based upon a given mammographic position. The 



blood volume and deoxygenation values are measured at a position closest to the 
mammographic position. On the other hand, in Fig. 15, the congruence of blood volume 
and oxygenation changes are used where ever it is observed in the breast. The use of the 
contralateral breast (i.e., the opposite breast) reduces the signals from the noncancerous 
5 tissue, but the measurement on the contralateral breast is not always feasible. 

Histograms shown in Figs. 16A, 16B, and 16 C plot the optical data shown in 
Figs. 14, 15 and 15A as a function of blood volume plus deoxygenation. 

Fig. 16A uses all the mammographic information for guidance and contains five 
false positives and zero false negatives. In Fig. 16B, I use the contralateral breast as a 
1 0 reference, but maintained congruence of the blood volume image and the deoxygenation 
image. Here, the four cancers form a significant cluster, again in the upper right hand 
portion with the data from patient 17 (i.e., metaplasia) included. Fig. 16B shows an 
improved segregation above 4 volts. The use of the contralateral breast clusters the 
normal patients below signals 3.5 volts, and gives a large gap between 3.5 V and 5.2 V. 
15 Furthermore, it clusters cancerous data above 5.2 V with one false positive and one false 
negative. Fig. 16C shows the optical data referenced to the breast model. However, I 
note that in the display using the breast model, many studies failed the congruence test 
(Fig. 16C); specifically, only 1 1 out of 30 met the criteria. This suggests that the use of 
the congruence test adds to the specificity. 
20 Quantitation of Co-registration. 

In principle, vasculogenesis (blood volume) and hypermetabolism (tissue 
hypoxia) occur in similar and often identical tissue locations. Therefore, I evaluated the 
congruence of the blood volume and deoxygenation images to further reinforce the 
identity of a suspicious mass. The following method for quantitation of the congruence 
25 was performed pixel by pixel. The first step was the normalization of the two images to 
equalize the maximum signals. A computer program was used to select the area and 
obtain the integrated value for the spatial congruence residual and for the blood volume 
signal. Then, subtraction pixel-by-pixel provided a residual value on which to base an 
estimate of the congruence of the two shapes, blood volume and deoxygenation. This 
30 was carried out for those shapes which appear by inspection to be congruent and the 
integral of the residual non-zero pixels is compared to the total signal. A simpler 



evaluation procedure can take the maximum value of the difference and divide it by the 

maximum value of the normalized value for the two images. 

The following test was applied to those images in the upper right region of Fig. 

14. Fig. 17 is a "four" dimensional graph summarizing 30 studies of patients with 

suspicious masses. The blood volume (Volts) was plotted on the abscissa and 

deoxygenation (Volts) on the ordinate. The measured size of image was depicted as the 

circle diameter and the percentage congruence between the blood volume image and the 

deoxygenation image is shown by a color scale. The following is the color coding of the 

percentage of congruence signals given by a color scale in Fig. 17. 

j _j maximum overlap residual 
{ maximum blood volume signal 

The diagram is summarized as follows: 

1 . The size of the image of a suspicious mass (plotted as one half its longest dimension). 

2. The congruence of blood volumes and blood deoxygenation plotted in color. 

3. The blood volume in the congruent region measured in Volts (scale of the abscissa). 

4. Blood deoxygenation in the congruent region (scale of the ordinate). 

For example, a red dot of a given size in the upper right hand quadrant of the 
diagram indicates a suspicious mass, having a large blood volume and a large 
deoxygenation with a high spatial congruence of the two parameters. Several 
measurements in this diagram either lacked congruence, lacked high blood volume, or 
lacked large blood deoxygenation. These are mainly in the lower left quadrant. 

This preliminary study has employed a combination of angiogenesis and 
hypermetabolism as indices of tumor regenesis. While the study has been guided by 
mammographic/ultrasound location of the tumor volume, only in one case has that 
information been required to identify the site at which the optical measurements were to 
be taken, the resolution of the images of blood volume and deoxygenation is adequate to 
determine their congruence. The congruence level is at the >70% level for the images of 
cancer. Thus, morphology is an asset of these studies and a 2D projection of the tumor 
size is indicated by the images. 

Localized increases of blood volume and of deoxygenation have been used as the 
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principle criteria of this study, others may add to the sensitivity/specificity. The first one 
of these is an estimate of congruence which is over 80% for those studies in the upper 
right hand corner of the diagram and poor for those in the lower left sector. Fig. 17 
indicates congruence of blood volume and deoxygenation with a color scale indicated on 

5 the margin of the figure. It is seen that the congruence is very high in cancer patients. 
Congruence is low elsewhere in the lower left corner and when it is less than 70% a blue 
color is used. The comparison of the non-congruent images is of interest, usually one 
feature or the other, blood volume or deoxygenation is well defined whilst the 
complementary feature is quite diffuse. A third parameter that can be added to the figure 

10 (and has been done in the case of the cancers) is half the size of the image in centimeters 
and most of the cancer bearing tissues areas are between 1 and 2 cm in diameter. 
Breast Matching Medium 

In other studies we have employed the breast immersed in a liquid of u a , fi s ' 
approximately matching the tumor, intending to render it optically invisible except for 

15 the inhomogeneity introduced by the tumor. This has been found to be quite important 
when transmission measurements are made through the breast. In the above study, 
however, each patient was lying on her back placing the breast against the chest, the 
body wall and ribcage, which created a matching background for the optical radiation. 
In the case of fatty breasts, the soft tissue "runs" over the chest and forms a thinner and 

20 optically more acceptable system, where the suspicious mass may be much better 

detected than in pendant breasts. The lymphatics could also be studied by this method. 
Higher Resolution Images 

In previous studies, optical tomography has attempted to mimic the X-ray image 
by a 2D projection of absorbance usually in 2 planes as shown, for example, by Fantini, 

25 et al. using a compressed breast and transmission time or frequency domain methods. 
The success of this techniques is based upon the ability of the radiologist to identify the 
structures of either scattering or absorbing material that differ from the normal breast. 
However, a high resolution is required to delineate such structural features on which 
mammographic identification of cancer is usually based, and high resolution is time 

30 intensive as well as apparatus intensive, i.e., numerous source detector combinations are 
required to achieve resolution comparable to PET/MRI. In this study, imaging resolution 
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is employed mainly to increase the signal to noise ratio in quantifying optical properties 
of the tumor with respect to the normal, often adipose tissue. 

In the above systems, imaging resolution is employed mainly to increase the 
signal to noise ratio in quantifying optical properties of the tumor with respect to normal 
5 tissue, or with trespect to a model of a normal tissue. However, the blood volume, 
oxygenation and deoxygenation data collected by the optical systems do not depend 
critically upon high resolution imaging. 

An optical system with an increased number of sources and detectors will render 
higher spatial resolution. Furthermore, a larger source-detector separation (i.e., the input 
1 0 port to detection port separation) achieves deeper penetration of the introduced optical 
radiation. By using selected separation values, the above-described imaging systems can 
collect three-dimensional optical data that are used for three dimensional reconstruction. 

Additional embodiments are within the following claims: 



- i» - 

Claims: 

1 . An optical system for in vivo, non-invasive examination of breast tissue of a 
female subject comprising: 

an optical module including an array of optical input ports and detection ports 
5 located in a selected geometrical pattern to provide a multiplicity of photon migration 
paths inside an examined region of the tissue, each said optical input port being 
constructed to introduce visible or infrared light emitted from a light source, each said 
optical detection port being constructed to receive photons of light that have migrated in 
the examined tissue region from at least one of said input ports and provide said received 
10 light to a light detector; 

a controller constructed and arranged to control operation of said light source and 
said light detector to detect light that has migrated over at least one of said photon 
migration paths; and 

a processor connected to receive signals from said detector and arranged to form 
15 at least two data sets, a first of said data sets representing blood volume in the examined 
tissue region and a second of said data sets representing blood oxygenation in the 
examined tissue region; said processor being arranged to correlate said first and second 
data sets to detect abnormal tissue in the examined tissue region. 

20 2. The optical system of claim 1 wherein said second data set includes 

hemoglobin deoxygenation values. 

3. The optical system of claim 1 wherein said processor is arranged to form 
a third data set being collected by irradiating a reference tissue region. 

25 

4. An optical system for in vivo, non-invasive examination of breast tissue 
of a female subject comprising: 

an optical module including an array of optical input ports and detection ports 
located in a selected geometrical pattern to provide a multiplicity of photon migration 
30 paths inside an examined region of the tissue, each said optical input port being 

constructed to introduce visible or infrared light emitted from a light source, each said 
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optical detection port being constructed to receive photons of light that have migrated in 
the tissue from at least one of said input ports and provide said received light to a light 
detector; 

a controller constructed and arranged to control operation of said light source and 
5 said light detector to detect light that has migrated over at least one of said photon 
migration paths; and 

a processor connected to receive signals from said detector and arranged to form 
at least two data sets, a first of said data sets being collected by irradiating an examined 
tissue region of interest and a second of said data sets being collected by irradiating a 
10 reference tissue region having similar light scattering and absorptive properties as the 
examined tissue region, said processor being arranged to correlate said first and second 
data sets to detect abnormal tissue in the examined tissue region. 

5. An optical system for in vivo, non-invasive examination of breast tissue 
15 of a female subject comprising: 

an optical module including an array of optical input ports and detection ports 
located in a selected geometrical pattern to provide a multiplicity of photon migration 
paths inside an examined region of the tissue or a model representing biological tissue, 
each said optical input port being constructed to introduce visible or infrared light 
20 emitted from a light source, each said optical detection port being constructed to receive 
photons of light that have migrated in the tissue or the model from at least one of said 
input ports and provide said received light to a light detector; 

a controller constructed and arranged to control operation of said light source and 
said light detector to detect light that has migrated over at least one of said photon 
25 migration paths; and 

a processor connected to receive signals from said detector and arranged to form 
at least two data sets of two tissue regions, a first of said data sets being collected by 
irradiating an examined tissue region and a. second of said data sets being collected by 
irradiating a region of a tissue model having selected light scattering and absorptive 
30 properties, said processor being arranged to correlate said first and second data sets to 
detect abnormal tissue in the examined tissue region. 
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6. The optical system of claim 1 or 3 further including a second optical 
module including an array of optical input ports and detection ports located in a selected 
geometrical pattern to provide a multiplicity of photon migration paths inside an 
examined region of the tissue, each said optical input port being constructed to introduce 
visible or infrared light emitted from a light source, each said optical detection port 
being constructed to receive photons of light that have migrated in the examined tissue 
region from at least one of said input ports and provide said received light to a light 
detector; said processor being arranged to receive optical data from both said optical 
modules. 

7. The optical system of claim 1, 4 or 5 wherein said processor is arranged 
to correlate said first and second data sets by determining congruence between data of 
said two sets. 

8. The optical system of claim 7 wherein said processor is programmed to 
order said first and second data sets as two-dimensional images and to determine said 
congruence using said two-dimensional images. 

9. The optical system of claim 7 wherein said processor is programmed to 
order said first and second data sets as two-dimensional images and to determine said 
congruence using the following formula: 



10. The optical system of claim 1, 2,3, 4, 5, 6, 7, 8 or 9 wherein said 
processor is further arranged to determine a location of said abnormal tissue within the 
examined tissue region. 




1 1 . The optical system of claim 1 , 4 or 5 wherein said processor is adapted to 
produce from said data set an image data set by implementing an optical tomography 



algorithm. 
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12. The optical system of claim 1 1 in which said optical tomography 
algorithm employs factors related to determined probability distribution of photons 
attributable to the scattering character of the tissue being imaged. 

13. The optical system of claim 1 , 4 or 5 wherein said controller is arranged 
to activate said source and said detector to obtain a first selected distance between said 
input and detection ports and said processor is arranged to form said data set for said 
first distance. 

14. The optical system of claim 13 wherein said processor produces an imge 
data set from said data set formed for said first distance. 

15. The optical system of claim 1 3 wherein said controller further is arranged 
to activate said source and said detector to obtain second distance between said input and 
detection ports and said processor is arranged to form another data set for said second 
distance. 

1 6. The optical system of claim 1 1 , 1 2, 1 3, 1 4 or 1 5 further including a 
display device constructed to receive said image data set from said processor and to 
display an image. 

17. The optical system of claim 1, 4 or 5 further comprising 

a first oscillator constructed to generate a first carrier waveform at a first 
frequency on the order of 10 8 Hz, said first frequency having a time characteristic 
compatible with the time delay of photon migration from said input port to said 
detection port; 

said light source being coupled to said first oscillator and constructed to generate 
said light modulated by said first carrier waveform; 

a phase detector constructed to determine change in waveform of the detected 
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light relative to the waveform of the introduced light and measure therefrom the phase 
shift of said detected light at said wavelength, said phase-shifted light being indicative of 
scattering or absorptive properties of the examined tissue region; and 

said processor being arranged to form said data set based on the measured phase 

shift. 

18. The optical system of claim 1 7 further comprising 

a second oscillator constructed to generate a second waveform at a second 
frequency; 

said detector being arranged to receive a reference waveform at a reference 
frequency offset by a frequency on the order of 10 3 Hz from said first frequency and to 
produce a signal, at said offset frequency, corresponding to said detected radiation; and 

said phase detector being adapted to compare, at said offset frequency, the 
detected radiation with the introduced radiation and to determine therefrom said phase 
shift. 

19. The optical system of claim 1, 4 or 5 further comprising: 

an oscillator constructed to generate a first carrier waveform of a selected 
frequency compatible with time delay of photon migration from said input port to said 
detection port; said light source being connected to receive from said oscillator said 
carrier waveform and constructed to generate optical radiation modulated at said 
frequency; 

a phase splitter connected to receive said carrier waveform from said oscillator 
and produce first and second reference phase signals of predefined substantially different 
phases; 

first and second double balanced mixers connected to receive from said phase 
splitter said first and second reference phase signals, respectively, and connected to 
receive from said detector said detector signal and to produce therefrom a in-phase 
output signal and a quadrature output signal, respectively; and 

said processor being connected to said double balanced mixers and arranged to 
receive said in-phase output signal and said quadrature output signal and form therefrom 



said data set. 
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20. The optical system of claim 19 wherein said processor is arranged to 
calculate a phase shift (0 X ) between said light introduced at said input port and said light 
detected at said detection port prior to forming said data set. 

21 . The optical system of claim 19 wherein said processor is arranged to 
calculate an average migration pathlength of photons scattered in the examined tissue 
between said optical input port and said optical detection port prior to forming said data 
set. 

22. The optical system of claim 21 wherein said processor further employs said 
pathlength in quantifying hemoglobin saturation (Y) of the examined tissue. 

23. The optical system of claim 19 wherein said processor is arranged to 
calculate a signal amplitude (A x ) determined as a square root of a sum of squares of said 
in-phase output signal and said quadrature output signal prior to forming said data set. 

24. The optical system of claim 23 further comprising: 

a narrow band detector connected to receive from said optical detector said 
detector signal and to produce a DC output signal therefrom; and 

said signal processor further determining a modulation index (M x ) as a ratio of 
values of said signal amplitude and said signal amplitude plus said DC output signal. 

25. The optical system of claim 1 , 4 or 5 further comprising: 

at least one oscillator constructed to generate a carrier waveform of a selected 
frequency, said light source being operatively connected to said oscilator 
constructed to generate light of a visible or infrared wavelength, said light being 
intensity modulated at said frequency to achieve a known light pattern; 

said controller constructed to control the emitted light intensity or phase 
relationship of patterns simultaneously introduced from multiple input ports, said 
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introduced patterns forming resulting radiation that possesses a substantial gradient of 
photon density in at least one direction, said resulting radiation being scattered and 
absorbed over said migration paths; 

said detector constructed and arranged to detect over time the resulting radiation 
that has migrated in the tissue to said detection port, and 

said processor being further arranged to process signals of said detected resulting 
radiation in relation to said introduced radiation to create said data sets indicative of 
influence of the examined tissue upon said substantial gradient of photon density of said 
resulting radiation. 

26. The optical system of claim 25 further comprising a phase detector 
constructed to detect the phase of the detected radiation and provide said phase to said 
processor. 

27. The optical system of claim 25 further comprising an amplitude detector 
constructed to detect the amplitude of the detected radiation and provide said amplitude 
to said processor. 

28. The optical system of claim 25 wherein the phase relationship of light 
patterns introduced from two input ports is 180 degrees. 

29. The optical system of claim 1 , 4 or 5 wherein said light source produces 
relatively long light pulses and the processor forms said data set by subtracting 
amplitude of two said pulses emitted from two input ports located symmetrically relative 
to one detection port. 

30. The optical system of claim 1, 4 or 5 constructed to introduce and detect 
photons at two wavelengths selected to provide sensitivity to a tissue constituent. 

3 1 . The optical system of claim 30 wherein said tissue constituent is an 
endogenous pigment. 



32. 

hemoglobin. 



-43 - 

The optical system of claim 3 1 wherein said endogenous pigment is 



33. The optical system of claim 30 wherein said tissue constituent in an 
exogenous pigment. 

34. The optical system of claim 33 in which said exogenous pigment is a 
selected contrast agent. 

35. An optical method for in vivo, non-invasive examination of breast tissue of a 
female subject comprising: 

providing an optical module including an array of optical input ports and 
detection ports located in a selected geometrical pattern to provide a multiplicity of 
photon migration paths inside an examined region of the tissue; 

placing said optical module on the breast of a female subject; 

introducing visible or infrared light from at least one said optical input port into 
an examined tissue region and receiving photons of light that have migrated in the 
examined tissue region to at least one of said detection ports; 

detecting said received photons by at least one optical detector optically coupled 
to said least one detection port; 

controlling said introducing and detecting steps to collect optical data 
corresponding to photons of light that have migrated between selected input and 
detection ports; 

processing said optical data to form at least two data sets, a first of said data sets 
representing blood volume in the examined tissue region and a second of said data sets 
representing blood oxygenation in the examined tissue region; and 

correlating said first and second data sets to detect abnormal tissue in the 
examined tissue region. 
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